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Three coupled Ginzburg-Landau equations for hexagonal patterns with broken chiral symmetry 
are investigated. They are relevant for the dynamics close to onset of rotating non-Boussinesq or 
surface-tension-driven convection. Steady and oscillatory, long- and short-wave instabilities of the 
Qv. ' hexagons are found. For the long-wave behavior coupled phase equations are derived. Numerical 

O^ . simulations of the Ginzburg-Landau equations indicate bistability between spatio-temporally chaotic 

patterns and stable steady hexagons. The chaotic state can, however, not be described properly with 
the Ginzburg-Landau equations. 
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O ' I. INTRODUCTION 
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r-^ , Convection has played a key role in the elucidation of the spatio-temporal dynamics arising in nonequilibrium 

^s^ • pattern forming systems. The interplay of well-controlled experiments with analytical and numerical theoretical work 

^\ ' has contributed to a better understanding of various mechanisms that can lead to complex behavior. From a theoretical 

^ "1 point of view the effect of rotation on roll convection has been particularly interesting because it can lead to spatio- 

O \ temporal chaos immediately above threshold where the small amplitude of the pattern allows a simplified treatment. 

I . Early work of Kiippers and Lortz MEJ showed that for sufficiently large rotation rate the roll pattern becomes unstable 

"{—{ ' to another set of rolls rotated with respect to the initial one. Due to isotropy the new set of rolls is also unstable and 

d , persistent dynamics are expected. Later Busse and Heikes Q confirmed experimentally the existence of this instability 

Mh and the persistent dynamics arising from it. They proposed an idealized model of three coupled amplitude equations 

J> , in which the instability leads to a heteroclinic cycle connecting three sets of rolls rotated by 120° with respect to each 

• ""j ■ other. Recently the Kiippers-Lortz instability and the ensuing dynamics have been subject to intensive research, both 

rS ' experimentally |Q-^ and theoretically p-[l5|| . It is found that in sufficiently large systems the switching between rolls 

^ of different orientation looses coherence and the pattern breaks up into patches in which the rolls change orientation at 

■ " " ' different times. The shape and size of the patches changes persistently due to the motion of the fronts separating them. 

Other interesting aspects induced by rotation are the modification of the dynamics of defects [Rm and an unexpected 

transition to square patterns ||l^ . 

In this paper we are interested in the effect of rotation on hexagonal rather than roll (stripe) patterns as they 
arise in systems with broken up-down symmetry (e.g. non-Boussines convection or surface-tension driven convection). 
Complex dynamics, if they are indeed induced by the rotation, are likely to differ qualitatively from those in roll 
patterns due to the difference in the symmetry of the pattern. Considering the small-amplitude regime close to onset, 
we use coupled Ginzburg-Landau equations. On this level the Coriolis force arising from rotation manifests itself as 
a breaking of the chiral symmetry. We therefore consider quite generally the effect of chiral symmetry breaking on 
weakly nonlinear hexagonal patterns. Since the equations are derived from the symmetries of the system we expect 
them to capture the generic behavior close to onset. 

The dynamics of strictly periodic hexagon patterns with broken chiral symmetry have been investigated in detail by 
Swift pql and Soward |jl3] . They found that the heteroclinic orbit of the Busse-Heikes model is replaced by a periodic 
orbit arising from a secondary Hopf bifurcation off the hexagons. Their results have been confirmed in numerical 
simulations of a Swift- Hohenberg- type model EQ]. The competition between hexagons, rolls, and squares in rotating 
Benard-Marangoni convection has been considered in pi| . In the present paper we focus on the impact of rotation 
on the side-band instabilities of steady hexagon patterns, i.e. on instabilities that introduce modes with wavelengths 
or orientation different than those of the hexagons themselves. Thus, we extend the work of Sushchik and Tsimring 
p2| to the case of broken chiral symmetry. We find that rotation can increase the wavenumber range over which the 



hexagons are stable with respect to long-wave perturbations. For larger values of the control parameter, however, 
additional short-wave instabilities arise. The long- and the short-wave modes can be steady or oscillatory. While in 
most cases they eventually lead to stable hexagon or roll patterns with different wavevectors, they can also induce 
persistent dynamics that can apparently not be described with Ginzburg-Landau equations. 

The paper is organized as follows. In the following section we use symmetry arguments to introduce the appropriate 
Ginzburg-Landau equations. The stability with respect to long-wave perturbations is addressed in section III in 
which the coupled phase equations for the system are derived. General perturbations (within the Ginzburg-Landau 
framework) are considered in section IV. In section V we investigate numerically the nonlinear behavior resulting from 
the side-band instabilities. Conclusions are given in section VI. 

II. AMPLITUDE EQUATIONS 

We consider small-amplitude hexagon patterns in systems with broken chiral symmetry. For strictly periodic patterns 
the amplitudes Ai of the three sets of rolls (stripes) that make up the hexagon satisfy then the equations pJ 



dtAi = eAi + aoA2A3 - giAi\Ai\^ - 92Ai\A2\'^ - 53^1^3 



(1) 



where the equations for the other two amplitudes are obtained by cyclic permutation of the indices and e is a small 
parameter related to the distance from threshold. The overbar represents complex conjugation. These equations can 
be obtained from the corresponding physical equations (e.g. Navier-Stokes) using a perturbative technique with the 
usual scalings {A ~ 0(e^/^) and dt ~ 0{e)). In order for all terms in (|l|) to be of the same order the coefficient of the 
quadratic term must be small, ao ~ 0(e^/^). This term arises from a resonance of the wavevectors of the three modes 
in the plane. The broken chiral symmetry manifests itself by the cross-coupling coefficients not being equal, 172 7^ 53- 

For completeness it should be noted that rotation leads in convection not only to a chiral symmetry breaking but 
also to a (weak) breaking of the translation symmetry due to the centrifugal force. In the following we will consider 
it to be negligible. In addition, for sufficiently small Prandtl number rotation can render the primary instability 
oscillatory |2|,|l|l. 

In order to analyze the possibility of modulational instabilities spatial derivatives must be included in Eq. (|l|) . We 
take the gradients in both directions to be of the same order, V ~ ©(e^^^) p5[, and retain both linear and quadratic 
gradient terms. After rescaling the amplitude, time, and space we arrive at the equations. 



dtAi = ^iAi + (m • V)Mi +A2A3 - Ai\Ai\^ -{iy + i>)Ai\A2\^ - {v ~ ^)Ai\A^ 
+i{ai + a)A2(n3 • V)1'3 + i(ai - a)A3(n2 • V)A2 
+ia2 (^2(1-3 • V)l3 - ^3(1-2 • V)32) 



(2) 



where now all the coefficients are 0(1), and n^ and t^ represent the unit vectors parallel and perpendicular to the 
wavenumber k^ (Fig. |l|). The cross-coupling coefficients have been rewritten in terms of v and v, with v being 
proportional to rotation and therefore giving a measure of the chiral symmetry breaking. In the gradient terms the 
chiral symmetry breaking manifests itself in the terms proportional to a. 



a) 




b) 



Kn 




FIG. 1. a) Critical wavenumbers ki,k2,k3. They satisfy the relation ki + k2 + ka = 0. Up to a global rotation '3/, the 
position of the wavevectors is given by the modulus of two of them and the angle between them. The symbols rii and n denote 
the unit vectors in the parallel and perpendicular directions to the wavevectors, respectively, b) The replacement Ai — > Aie*^"' 
represents a change in the magnitude of the wavevector while Ai -^ A^e*^^' represents a rotation of it. 

The influence of the nonlinear gradient terms in (0) involving ai and a2 has been studied by several authors 1 26-281. 
The origin of the new term involving a is best understood by considering the coefhcient a of the quadratic term. The 
gradient terms arise from its dependence on the wavenumber of the modes involved, a = a(ki,k2,k3) i.e. when the 
equation is considered in Fourier space. Due to the resonance condition (ki +k2 +k3 = 0) we can drop the dependence 
on one of the wavevectors. Then, up to an arbitrary global rotation $, the system can be specified by the angle Q 
between the wavevectors k2 and k3 and their moduli ^2 and k^ (see Fig. |l|). Due to isotropy the coefficient a cannot 
depend on the global rotation and can therefore be expressed as a = a(fc2,fc3,6). When evaluated at the critical 
values of the wavenumbers ao is given by ao = a{k2, k^, Q'^ = 27r/3) (in Eq. (0) we take the normalization ag = !)■ 
Since a change in the modulus of k^ can be effected by the replacement Ai -^ Aie^^'^' the dependence of a on the 
moduli fc2 and k^ can be represented in real space by 

dct —r , „s—r da —r , „s—r ,„S 

--^3 n2-V A2 + T^A2 ng-V ^3- (3 

ok2 oks 

When the chiral symmetry is broken da/dk2 and da/dk^ need not be equal and it is convenient to introduce the 
coefficients 

1 / da da\ ^ _ 1 f da da\ ,^. 

^' = 2[dk-, + dk-2) ^''^'' = 2[wrdk-2) ^'^ 

with ai even and a odd in the amplitude of the symmetry breaking. 

On the other hand, a variation in the angle between k2 and k3 is represented by 

^^^(A2{^3-V)A,-Mt2-V)A2) (5) 

with only one coefficient a2 — {da/dO)/kc. This term is invariant under reflections interchanging modes A2 and 
A3, since in contrast to the normal vector n^ the tangential vector Ti changes sign in this reflection. Therefore the 
coefficient a2 is even in the amplitude of the symmetry breaking. 

Equation (0) admits hexagon solutions Ai = jJe^^'^i-^+^'t'i vifith a sligthly offcritical wavenumber (ki = k^ + K;, 
K <^ kc), with 



H= ^^^^ . , $ = <?;'i + 02 + 03 = 0. (6) 



The stability of this solution to perturbations with the same wavevectors has been studied by several authors |18yi9(l 
and can be summarized in the bifurcation diagram shown in Fig. 0. The hexagons appear through a saddle-node 
bifurcation at /i = /is„, 

^-^' 4(1 + 2.) +^' ^^^ 

and become unstable via a Hopf bifurcation at fi — fin , 

^" = {u-iy + ' ^^> 

with a critical frequency ujc — 2^/?>i>{l + 2Kai)^ /{v — 1)^. Note that the Hopf frequency does not depend on a. The 
Hopf bifurcation is supercritical and for /i > ^h stable oscillations in the three amplitudes of the hexagonal pattern 
arise with a phase shift of 27r/3 between them [|l8|-p0|i resulting in what we are going to call oscillating hexagons. As 
IX is increased further, eventually a point 11 — iihet is reached at which the branch of oscillating hexagons ends on the 



branch corresponding to the mixed-mode solution in a global bifurcation involving a heteroclinic connection. Above 
this point the only stable solution is the roll solution whose stability region is bounded below by 



1 



fJ-R 



{v + 1/ - l){v - V -1) 



(9) 



When \v\> V — 1 the rolls are never stable and the limit cycle persists for arbitrary large values of ji. In the absence 
of the quadratic terms in Eq. ffl) this condition corresponds to the Kiippers-Lortz instability of rolls. 
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FIG. 2. Sketch of a bifurcation diagram for a fixed value of the wavenumber K. 



III. LONG- WAVE APPROXIMATION: PHASE EQUATION 



Already below the Hopf bifurcation to oscillating hexagons the hexagonal pattern can be unstable to side-band 
perturbations. The behavior of long-wavelength modulations is described by the dynamics of the phase of the periodic 
structure. The phases (pi of the three modes of the hexagonal pattern can be combined to define a phase vector 
4> = {4'x,4'y): the components 0j; = —(02 + ^s) and (j)y = (02 — 4'3)/V^ of which are related to the translation modes 
in the x - and the y -direction. In the chirally symmetric case the phase vector satisfies the coupled diffusion equations 
1290 



Todtct} = D^_W^4> + {D« - D^)W{W ■ 0), 



(10) 



and can be decomposed into a longitudinal (irrotational) and a transversal (divergence free) part (p = Vipi -f V x ez^t- 
The fields V'/,t each satisfy a diffusion equation with diffusion constants Dy and D_l, respectively. 

We can use symmetry arguments to derive the form of the phase equation when the chiral symmetry is broken. We 
consider a general diffusion coefficient which is a tensor of rank four, 



df 



_ £)jkl 



djd, 



'k(pi 



(11) 



with summation over repeated indices implied. Invariance under reflection and rotations of 60° restrict the number 
of possible independent coefficients Df . In the case of broken chiral symmetry we split the diffusion tensor into two 
parts: one even under reflections, the other odd. 



jjiki^-^,ki 



nTDi 



kl 



where il™ is the antisymmetric tensor of rank two given by 



Q = 



-u. 
oj 



(12) 



(13) 



with uj giving the strength of the chiral symmetry breaking (e.g. the rotation frequency). D and D are even functions 
of UJ. As generators of the symmetry group we can take rotations of 60° and reflections in x, 



Requiring that D and _D are invariant under the operations (|l4| ) one can show that the most general form of the phase 
equation with broken quiral symmetry is given by 

dtct> - Z?^V2</,+ (Z?|| - Z?^)V(V • 0) - i?xi(e. X V20) + i?x.(e. x V)(V • 0), (15) 

where e^ is a unit vector in the direction perpendicular to the plane. 

It is worth emphasizing that, although the coefficients of this equation can be derived from the amplitude equations, 
its form is given by symmetry arguments and is, therefore, generic and valid even far from threshold. To derive the 
phase equation from the amplitude equations (0) we consider a perfect hexagonal pattern with a wavenumber slightly 
different from critical {k = kc + K) and perturb it, both in amplitude and phase, Ai = {H + T-^^e*^"''^"'"*'^'. Away 
from threshold, from the saddle-node and the Hopf bifurcation, the amplitude modes r'i,r2,r3, and the global phase 
$ = (/)i + 02 + 03 are strongly damped and can be eliminated adiabatically. Following the usual procedures (e.g. |3l|| ) 
we arrive at Eq. (nsh with 

D^ = \+ ^^ A \H^u\{ai + V?,a2f + -ia^]-\/?,HujaK-uKA, (16) 

4 u^ + w^ [_ 4 J 

D«=D^ + ]--- (n^aiiai - V3a2) - H{3ai - V3a2)K + 2K^] , (17) 

" 2 V I J 

Dxi = A A \u;H^[iai + V3a2f + M^] + V3Hu&K-ujKH, (18) 

M^ + ti>^ 1^ 4 J 

Dy,^^-(V3H'^ai-V3HK\. (19) 

where 

to = 2VSH^i>, (20) 

u = 2H'^{l-iy) + 2{l + 2Kai)H, (21) 

V ^2H^{l + 2i')-{l + 2Kai)H. (22) 

The coefficients Z?xi and D^^ are odd in the symmetry-breaking terms v and a. At the Hopf bifurcation curve u = 
implying H ~ {I + 2Kai)/{v — 1) and u) — lOc, while u = represents the saddle-node instability. 
Expanding the phase in normal modes <p — (p e**^''+'^* we obtain the dispersion relation 

a^ + (D\\ + D^)Q^a + {D\\D^ + Dy,,{Dy,, + i?x.)) Q^ = 0, (23) 

whose eigenvalues are 

1 r. 



'^M--^ 



D.+D^±J{D.-Di_Y-AD^,{D^,+D^,,) Q'. (24) 



When Dxi = the eigenvalues become simply ai — —D\\Q^ and (T2 = —D±Q^, corresponding to the eigenvalues of 
the irrotational and the divergence-free phase modes, respectively. If the rotation rate is small {Dx^, 0^2 "^ D±, Dy) 
we can expand ( |24| ) and obtain. 

For Dxi,Dx2 ^ Dn,D± this approximation is not valid and the longitudinal and transverse perturbations become 
coupled. An important novelty in this case is that the phase instability can become oscillatory. It occurs when the 
following conditions are satisfied, 

D\\+D±=0, (26) 

iDn-D^f-ADx,{Dx,+Dx2)^0- (27) 



In Fig. ra and Fig. (below) we represent the phase instabihty curves for a number of cases. For smaU values of 
the rotation rate zJ, the phase stability diagram is similar to that obtained in the absence of rotation, especially for 
small fx. As \K\ is increased, both real eigenvalues in (p3) go through zero consecutively as indicated by the dashed 
and solid lines in Fig.g. As /i is increased towards the Hopf bifurcation (/i -^ jih) the two lines merge and the phase 
instability becomes oscillatory as indicated by the solid lines. Note that, in contrast to the chirally symmetric case, 
the left and right stability limits do not merge at iiT = as the transition to oscillating hexagons is reached. Instead, 
they are open and over a range of wavenumbers the hexagons remain stable with respect to long-wave perturbations 
all the way to the Hopf bifurcation at fi — fin- Furthermore, while in the chirally symmetric case the analog of the 
Hopf bifurcation is transcritical (and steady) and leads discontinously to rolls, the Hopf bifurcation is supercritical 
and leads to oscillating hexagons M-EQ] . 



6.0 



1^ 



0.0 




Oscillating Hexagons 



i? Stable 
' 'f Hexagons ,\: ^ 




1^ 



0.0 



b) 



Oscillating Hexagons 




-1.0 -0.5 0.0 

Wavenumber K= 



0.5 

k-k 



-0.5 0.0 0.5 

Wavenumber K=k-k. 




-0.5 0.0 

Wavenumber K= 



FIG. 3. Stability diagrams for u — 2, a — ai — a2 = and a) !> = 0.05, b) !> = 0.2 and c) i> — 0.5. The dotted and 
short-dashed lines indicate the saddle-node and the Hopf bifurcation to oscillating hexagons, respectively. The solid and dashed 
lines correspond to the long-wave instabilities, while the open circles are the results from the general stability analysis. Rolls 
are stable with respect to hexagons above the dot-dashed line. 



IV. GENERAL STABILITY ANALYSIS 



We now consider arbitrary perturbations of the hexagonal pattern Ai = {H + a.ie**^''+'^*)e'^"*'', with ai, 02, as 
complex and solve the resulting 6x6 linearized system. Two of the six eigenvalues correspond to the global phase $ 
and the overall amplitude involved in the saddle-node bifurcation. In the regime of interest both are strongly negative. 
The next two correspond to the translation modes, and can be real or complex. It turns out that these modes can 
destabilize the hexagons not only via the longwave instabilities (p3) but also via short-wave instabilities as illustrated in 
Fig. Ha, where the solid and dashed lines correspond to the real parts of the complex and real eigenvalues, respectively. 
Finally, there is a pair of complex conjugate eigenvalues corresponding to the Hopf bifurcation to oscillating hexagons. 
For some parameter values these eigenvalues merge with the ones corresponding to the phase modes and their real 
parts become positive (Fig. 0b and 0c). Note that in the chirally symmetric case the eigenvalues are always real and 
the instabilities long- wave p2|. 
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FIG. 4. Growth rates (with i/ = 2, D = 0.1, ai = a2 = 0) when a) the translation modes are short-wave unstable (a = 
and /i = 2.54, q = 0.53). The Hopf modes interact with the translation modes and their eigenvalues become positive at finite 
wavenumber {a — 0.4): b) Just below the branch switching {fi = 3.49, q — —0.6), and c) when the Hopf mode is unstable 
(/i = 3.75, q = —0.57). The solid and dashed lines correspond to complex and real eigenvalues. The points at which these 
dispersion relations are obtained are marked by an asterisk in Figs. W and BI 

In what follows we will consider ai = a2 = for simplicity. Although non-zero values for these coefficients change 
the stability boundaries quantitatively, they are not found to induce any qualitatively different instability. 

In Figs. and a we present the stability limits obtained for several values of i) and a. The short-dashed and the 
dotted lines correspond to the Hopf and saddle-node bifurcations, respectively, while the dot-dashed line is the curve 
above which the rolls become stable with respect to the hexagons. We do not address their side-band instabilities. 
The circles correspond to the results of the general stability analysis while the solid and dashed lines are the stability 
limits in the long- wave approximation, as given by (p4h (with eigenvalues either real or a complex conjugate pair). 
The solid circles in Fig. p^ correspond to instabilities at finite wavenumber due to the Hopf modes. 

Fig. shows the stability limits for a = but C' ^ 0, i.e. the chiral symmetry is only broken at cubic order. While 
for small values of the control parameter the long-wave analysis gives the correct stability limits, for larger /z a steady 
short-wave instability preempts the long-wave instability before it becomes oscillatory. In this case the eigenvalues are 
complex for Q ^ 0, but they split into two real eigenvalues for larger Q, one of which becomes positive (Fig. Ek). As 
the coefficient v is increased the region in which the long-wave instability is relevant decreases (Fig. ^p) and shrinks 
to almost zero (Fig. ^c). 

When the quadratic gradient terms are different from zero the instability regions become asymmetric with respect 
to K = 0. The system is, however, invariant under the change a — > —a, K —^ —K, a^ —f —ai, and we will therefore 
consider only positive values of a. For small a the region of the steady long-wave instability becomes smaller at one 
side of the bandcenter and larger at the other (Fig. @a) . As a is increased, this region shrinks to zero for K < and 
disappears above the Hopf curve for i^ > 0. At this point all the instabilities are oscillatory (Fig. ||c). For a — 0.4 
and a = 0.7 (Fig. &,c) the stability limit is entirely given by the longwave results for K > 0, but for still larger values 
of a it becomes shortwavelength (Fig. |p). For K < there is a large region in which the instability is short-wave 
and oscillatory. Close to the Hopf bifurcation the instability switches from the translation modes to the Hopf mode 
(cf. Fig. m. As a is increased the region in which the instability is due to the Hopf modes grows. 
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FIG. 5. Stability diagrams for u = 2, i> — 0.1, ai — a2 — and a) a = 0.1, b) a = 0.4, c) a — 0.7 and d) a = 1.5. The 
solid circles correspond to the instabilities at finite wavenumber due to the Hopf modes. The other lines are as in Fig. H. The 
diamond in d) corresponds to the simulation in Fig. pi 



V. NUMERICAL SIMULATIONS 



In order to study the nonlinear behavior arising from the instabilities, we have performed numerical simulations of 
Eqs. (y). A Runge-Kutta method with an integrating factor that computes the linear derivative terms exactly has 
been used. Derivatives were computed in Fourier space, using a two-dimensional fast Fourier transform (FFT). The 
numerical simulations were done in a rectangular box of aspect ratio 2/v3 with periodic boundary conditions. This 
aspect ratio was used to allow for regular hexagonal patterns. 

We start with a perfect hexagonal pattern with a wavenumber in the unstable region and add noise. In all the 
cases we have considered the numerical simulations reproduce correctly the linear stability limits. Over most of the 
parameter regime the nonlinear evolution of the instabilities is qualitatively very similar. The perturbation growths 
(with or without oscillations, depending on the kind of instability) until it destroys the original hexagon pattern and 
then settles down to a stable periodic pattern. Therefore all the instabilities appear to be subcritical. Furthermore, it 
seems to be irrelevant whether the instability comes from the translation or the Hopf modes. The branch switching 
does not change the value of the unstable wavenumber nor the frequency (cf. Fig. Qb,c) and the behavior of the 
dispersion relation at lower values of the perturbation wavenumber does not play a role. For values of the control 
parameter for which rolls are unstable, the instabilities lead to a rotation of the original hexagonal pattern and to a 
change in its wavelength. Usually penta-hepta defects appear in the process. In the presence of rotation they annihilate 
each other quite fast, yielding a perfect pattern as the final state. For larger control parameters rolls become stable 
and the side-band instabilities of the hexagons eventually lead to roll patterns, independently of the specific type of 
the instability. 

For certain parameter values, however, more complicated behavior is found. This is shown in Fig. |6|, where we 
represent a reconstruction of the hexagonal pattern, '^ = J2i=i^i^^^'^ (top panel), as well as the corresponding 
Fourier spectrum of the amplitude Ai, Ai{K) (bottompanel). In this case the instability develops close to the initial 
wavenumber of the unstable hexagonal pattern (Fig. ^a). As time progresses, however, modes with ever increasing 
y-component of the wavevector are excited. Independent of the maximal wavevectors retained in the simulations 
(15.3 < Kmax < 46 with system size in the range 17.5 < L < 52.5) eventually the wavevectors with the largest 
possible y-components are excited and the peak in the spectrum displayed in Fig. ^ reaches the top border of 
the figure. Then the peak reemerges at the bottom border again, i.e. the wavevectors have very large negative y- 
components. This is shown more clearly in Fig. R where a cross-section of the Fourier spectrum in the y-direction for 
Kx = is shown for three times. At time t ~ 90 most of the excited modes have already reemerged at (large) negative 
values of Ky. Obviously, in these simulations the solutions cease to be numerically resolved already well before t — 60. 
Within the Ginzburg-Landau equations (0) the curve of marginal modes corresponds to a vertical line in the Fourier 
spectrum in Fig. M. Thus, the excited modes lie predominantly along the critical curve and the numerically observed 
behavior suggests that the correct evolution of the pattern would involve a trend towards a rotation of the pattern 
and a spreading of the Fourier modes over the circle of marginal modes. Such dynamics reflect explicitly the isotropy 



of the system. They cannot be captured within the Ginzburg-Landau equations, which break the isotropy through the 
choice of the wavevectors corresponding to the ampHtudes Ai. To represent dynamics as suggested in Fig. |6| correctly, 
models that retain the isotropy have to be used. This motivates the use of Swift-Hohenberg-type models |32p|, |20|]P| . 
Investigations of the complex dynamics that can arise from instabilities identified here have been performed in [34|. 
They show indeed a bistability between the ordered hexagons and a spatio-temporally chaotic state with an almost 
isotropic Fourier spectrum. 




FIG. 6. Numerical simulations of Eqs. (El) in a box of length L — 35 {Kmin = 0.1795) with 128 x 128 Fourier modes, for 
the values of the parameters v = 2, v = 0.1, ai = 02 = 0, a = 1.5, ^ = 3 and K = —4,Kmin (cf. Fig. Hd). The top panels 
show a reconstruction of the hexagonal pattern ^ = X/i=i ^is'''' " with kc = 2QKmin for the times t = 32, t = 52 and t = 56. 
In the bottom panels the corresponding Fourier spectra of the amplitude A\ are shown. Note the shifting of the peak of the 
distribution to ever larger y-components of the wavevector. 
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^Formally, isotropy can be recovered by modifying the gradient terms in Eq. 
amplitudes be allowed to vary rapidly in space. 



However, this requires that the 



FIG. 7. A cut of the Fourier spectrum (in logarithmic scale) of the amplitude Ai at K^ = for the same parameters as in 
Fig. o. As time progresses the peak of the distribution evolves towards higher wavenumbers and eventually reemerges at the 
negative end of the spectrum. 

VI. CONCLUSION 

In this article we have analyzed the effect of chiral symmetry-breaking on the stability of hexagonal patterns. Such 
patterns arise, for instance, in non-Boussinesq Rayleigh-Benard convection and in Marangoni convection, where the 
chiral symmetry can be broken by rotating the system. Focussing on the regime near threshold we have used the 
appropriate Ginzburg-Landau equations for the three modes making up the hexagon pattern. The chiral symmetry 
breaking introduces an asymmetry between the cubic coupling coefficients as well as a new nonlinear gradient term. The 
general linear stability analysis of these equations revealed long-wave as well as short-wave instabilities. The long-wave 
instabilities, which are captured with coupled phase equations, can be steady or oscillatory. For all parameter regimes 
investigated, the short-wave instabilities arise for larger values of the control parameter, but below the transition 
to oscillating hexagons. They can be due to the translation or the Hopf modes. In the latter case they are always 
oscillatory. 

In contrast to the Kiippers-Lortz instability of stripe patterns p] , no regime was identified in which hexagon patterns 
become unstable at all wavelengths. Nevertheless, persistent irregular dynamics of disordered hexagon patterns can 
apparently arise from the short-wave instability. Our numerical simulations of the Ginzburg-Landau equations indicate 
that the nonlinear evolution ensuing from the instability tends to introduce modes with wavevectors covering the whole 
critical circle. Of course, such a state in which the Fourier modes are distributed almost isotropically over the critical 
circle cannot be described by Ginzburg-Landau equations, since they break the isotropy at the very outset. This 
suggests the use of Swift-Hohenberg-type equations, which preserve the isotropy of the system. They are often used 
as truncated model equations (e.g. |3q|) to study the qualitative behavior of systems, but can under certain conditions 
also be derived from the basic (fluid) equations as a long- wave description p6|~^ . Recently, in such investigations of 
hexagons with broken chiral symm etry spatio-temporally chaotic states have been found to arise from the corresponding 
oscillatory short-wave instability ||34| . As in our simulations of the Ginzburg-Landau equations the spatio-temporal 
chaos persists although for the same parameters there exist also stable ordered hexagon patterns. This bistability 
is somewhat reminiscence of the coexistence of spiral-defect chaos and ordered roll convection in Rayleigh-Benard 
convection without rotation |39] . In the Swift-Hohenberg model the oscillatory short-wave instability can also lead to 
a supercritical bifurcation to hexagons that are modulated periodically in space and time [Q . No such state could be 
identified in the Ginzburg-Landau equations discussed here. 

Our results suggest that rotation may induce irregular dynamics in hexagonal convection patterns quite close to 
threshold. So far, disordered hexagon patterns (without broken chiral symmetry) have been found in Marangoni 
convection pt3[ far from threshold and also in experiments on chemical Turing patterns [ pd[ . In the latter case they 
appear to be due to the competition with the stripe pattern in a bistable regime. 

From previous work it is well known that the chiral symmetry breaking delays the transition from hexagons to 
stripe patterns. More specifically, the steady bifurcation to the unstable mixed state is replaced by a Hopf bifurcation 
to a state of coherently oscillating hexagons [|l8| 20[ . Their side-band instabilities can be investigated with the same 
Ginzburg-Landau equations as discussed here 
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